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Abstract 

A completely new approach to the Ising model in 1 to 5 dimensions is developed. We employ 
p, g-binomial coefficients, a generalisation of the binomial coefficients, to describe the magnetisation 
distributions of the Ising model. For the complete graph this distribution corresponds exactly to the 
limit case p = q. We take our investigation to the simple d-dimensional lattices for d = 1, 2, 3, 4, 5 
and fit p, g-binomial distributions to our data, some of which are exact but most are sampled. 
For d = 1 and d = 5 the magnetisation distributions are remarkably well-fitted by p, g-binomial 
distributions. For d = 4 we are only slightly less successful, while for d = 2, 3 we see some deviations 
(with exceptions!) between the p, q-binomial and the Ising distribution. We begin the paper by 
giving results on the behaviour of the p, g-distribution and its moment growth exponents given a 
certain parameterization of p, q. Since the moment exponents are known for the Ising model (or 
at least approximately for d = 3) we can predict how p, q should behave and compare this to our 
measured p, q. The results speak in favour of the p, g-binomial distribution's correctness regarding 
their general behaviour in comparison to the Ising model. The full extent to which they correctly 
model the Ising distribution is not settled though. 
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I. INTRODUCTION 



Choose a graph, e.g. a square lattice, on n vertices and compute its Ising partition 
function Z, keeping track of its terms according to their magnetisation, so that Z = Z0+Z1 + 
■ ■ ■ + Zn- The quotient Z^/Z is the probabihty of having k negative spins, or magnetisation 
M = n — 2k. The controlhng parameter of the partition function is the temperature, 
that is, for any given temperature the partition function provides us with a distribution of 
magnetisations. At infinite temperature this is simply the binomial distribution. At zero 
temperature, on the other hand, we receive a distribution with two peaks, one at k = 
and one at k = n, both with 50% of the probability mass. The distribution is always 
symmetrical. What happens between these two extreme temperatures? Though there are 
exceptions to the rule, for most graphs the distribution begins its life at high temperatures as 
a unimodal distribution with the peak at the middle k = n/2. As we lower the temperature 
the distribution gets increasingly wider until we reach a temperature where the distribution 
changes from unimodal to bimodal. Near this temperature, slightly above and slightly below, 
the distribution is particularly wide. Lowering the temperature even further the distribution 
develops two sharp peaks, both essentially gaussian, and the peaks move outwards. 

This article tries to model the distributions using p, g-binomial coefficients. In one case, 
the complete graph, they model exactly the distributions and for rf- dimensional lattices the 
dimension seems to determine how well they fit to the Ising distributions. Especially in the 
case ofd = 1 and d = 5 the Ising distributions are particularly well-fitted by the p, g-binomial 
coefficients. 

The paper begins in section [III by providing the necessary basic tools, such as the 
Pochhammer and g-Pochhammer symbol, g-binomial coefficients and finally the p, g-binomial 
coefficients. Some nice results on their properties are also stated, even though they are of 
no direct use to us in the rest of the paper. They are merely intended to give the reader 
a feel of how p, g-binomials behave. In section IIIII we define the p, g-binomial distribution 
and give an algorithm for finding values of p and g when the distribution is given as input. 
The problem is to determine an optimal choice for p and g. As it turns out we only have 
to focus on the value of the probability and the location of the distribution's peaks, at least 
for a bimodal distribution. The unimodal distribution always has its peak at the middle so 
in this case we instead take the quotient between the two middle probabilities as controlling 
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parameter. This quotient is unfortunately rather sensitive to noise, making it difficult to 
determine the parameters p and q for sampled data. 

Section [IV] gives detailed results for the special case when p = q. This case corresponds 
exactly to the complete graph and is the only case where we can give asymptotically exact 
expressions for the sum of the coefficients. In section |V] we provide some useful tools for 
working with p, g-binomial coefficients in the case when p ^ q. After this build-up of tools 
we are, at long last, ready to give some general results on the distribution of p, g-binomial 
coefficients in section IVTl Using the parameterization p = 1 + y/n and q = 1 + z/n we find 
the asymptotic value of y, given z, where the distribution is fiat in the middle. We also allow 
for a small change in y, using a higher order parameter a, so that we can follow properly how 
the distribution changes from unimodal to bimodal. However, the computations that we rely 
on involve some rather complicated series expansions that were made using Mathematica. 
These are much too long to fit into this paper. We have prepared a simplified Mathematica 
notebook that performs all the necessary computations. The interested reader can obtain it 
by contacting the ffist author. 

Section IVIII looks into the case of moving y and z along a line with any given slope as 
opposed to the previous section where z stays fixed. In section IVlllI we give exact scaling 
formulae for the moments of the distributions depending on the parameters a and z. For a 
given moment of these distributions we always obtain the same exponent on n, regardless 
of a and z. In section [IX] we try to remedy this by letting the previously fixed parameter z 
depend ever so slightly (at most logarithmically) on n. This is based on the assumption that 
the previous formulae in section IVTl still hold. However, we can now adjust the exponent of 
n though this comes at the cost of an extremely slow convergence. 

Section |X] defines the Ising model, laying the ground for studying distributions of mag- 
netisations, the intended application of our endeavour. In section |XT] we apply our tools 
to the d- dimensional lattice graphs for d = 1,2, 3, 4, 5 fitting p, g-distributions to simulated 
distributions and comparing them. 

A condensed reading, more suitable to the reader who is pressed for time, should include 
a look at (P, ([II]), ([271), ([2H]), (122]), (in]), (I^ necessary definitions and results 

concerning the basics. After that, the most important results are stated in equations (1621) . 
fl^S]) . fl^Bl) . (I109p and (I112p . After looking up the basic definitions regarding the Ising model 
the reader can skip to (11211) . In Section jXIj the reader can now pick and choose his favourite 
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lattice and look at the pictures. 



II. DEFINITIONS, NOTATIONS, THE VERY BASICS 



The g'-binomial coefficient 



\ — qn-i+l 



, qj^l,0<k<n 



-' 1 i=l 

is a natural extension of the standard binomial coefficient 



n 



ni 



-, < A; < n 



^kj k\{n-k)V 
The Pochhammer symbol (or shifted factorial) is defined as 

n— 1 



i=0 



SO that (1)„ = n\. Its g-deformed relative, the gr-Pochhammer symbol, is defined as 



n-l 



i=0 

The g-binomial coefficient can then be expressed as 

k\g {q; q)k (q-, Q)n-k g)fc 

The g-numbers are defined for any real number a as 



1- q"- 



and it is easy to show that 



limfalo = a 



Note that for integers n > 1 we have 



l + gH q^l 



1-q 



so that [n\q ^ n when q ^ 1. To continue, the g- number factorials are defined as 

n 



and then obviously 



limfnl 

9^1 



The g-binomial can now be defined in an alternative way as 

\n\' 



and then it follows that 



lim 



n 






_k_ 




(:) 



Quite analogously it is easy to verify that 



lim 



(10) 



(11) 



(12) 



(13) 



Finally, note also that 



9^1 (1 - g)" ' ' 
can be viewed as a formal polynomial in q of degree k{n — k) 



where the coefficient of counts the number of /c-subsets of {1, . . . ,n} with element sum 
j + k {k + l)/2. It is thus a polynomial with positive coefficients. 

We have so far only stated what belongs to the standard repertoire on the subject. The 
g-binomials and the g-Pochhammer function have many interesting properties and we point 
the interested reader to the books^,-,- and especially the charming little book^. For more on 
the standard binomial coefficient and Pochhammer function we recommend^ which contains 
a wealth of useful information. 

A natural extension of the g-binomial coefficient, the p, g-binomial coefficient, was defined 



in^ as 



n 



p 



n—i+l 



- q 



n—i+l 



V 1^ 1) < k < n 



(14) 



n 




n 




_k_ 




_k_ 





p' — g' 

Clearly, in the case p = 1 this reduces to a g-binomial coefficient. Also, note that p and g 
are interchangable so that 

(15) 



(16) 



It is an easy exercise to show the following identity and we leave this to the reader. 



Just as the standard binomial coefficients, their p, g-analogues are also symmetric 



n 




n 




k_ 


p.? 


n — k 





n 
k 



p 



k {n—k) 



P>9 



q 



k {n~k) 



J q/p 



(17) 



J p/q 



As a corollary it follows that 



lim 

p,g-s>r 



„fc(n-fc) 



(18) 



Any identity involving g-binomial coefficients can then be extended to a p, g-binomial 
identity by ffist replacing q with q/p and then use the identity f|T7|l . For example, the 
binomial theorem 

k=0 

has the g-analogue 



X 



n-l 



n(i + -<;')=E 

£=0 k=0 

Now replace q with q/p and then use the identity 

— k in—k) 



q\2j X 



P 



q/p 



from f ll7p above. This gives us 



n-l 



n (1 + X (q/pY) = J2 



1=0 



p--k(n^k) ^qjp^iX) 



k=o -'p^q 

After multiplying both sides with ^(2) this simplifies into 



n-l 



f n—k\ / k\ 



p\ 2 J q\2j X 



p,q 



n(p'+-9')=E 

£=0 k=0 

which was also shown in^ using a recursion technique. 

As another application we consider the Chu-Vandermonde identity 



m + n 
k 



£=0 



m \ in 

k-l 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



which is a fairly direct consequence of the binomial theorem applied to the product (1 + 
x)"^ (1 + x)". The g-Vandermonde identity, see e.g.- and (though misprinted)^, can be stated 

as 



m + n 
k 



E 



-I q £=0 



k 



n 



£{m-k+£) 



(25) 



Using ( !T7|) above we can now obtain a p, g-analog of this. 
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Theorem II. 1. 



m + n 
k 



m 
k- 



n 



(26) 



Proof. In the g-Vandermonde identity (l25l) . replace q with q/p and muhiply both sides with 
pk{m+n-k)^ Using f|T7|l the left hand side is now a pure p, g-binomial coefficient. The ith. 
term of the right hand side is 



m 




n 


_k-i_ 


q/p 


i_ 



p 



k{m+n-k) (^g ^ pY (^~^+^) 



<i/p 



m 




n 


_k-i_ 


q/p 





pk {m+n-k)-e (m-k+e) (m-k+e) 



q/p 



P 



(k-e) (m-(k-e)) 



m 



k 



P 



e (n-e) 



m 
k- 



-I P,Q 



q/p 
n 



n 



pin-£) (k-e) (m-k+e) 



q/p 



p{n-e) (k-e) ^e (m-k+e) 



and the theorem follows. 



□ 



The g-binomial coefficients have been shown to form a log-concave (and thus unimodal) 
sequence for q > 0, see e.g.- and^. However, for the p, g-binomial coefficients this does not 



always hold. Rewriting them as a product like in ( 1T7|) we have in fact a product of two 

n 

for k = 0, . . . ,n. The ffist sequence is log-concave 

q/p 



sequences; that of p'^*^" and 



k 



for p > 1 and log-convex for p < 1. It is well-known that the element- wise product of two 
log-concave positive sequences is also log-concave. So, if p > 1 and g > then the sequence 
of is log-concave. 

>J p,q 

We conjecture that for p,q > the sequence can be either unimodal, with the maximum 
ai k = \n/2\, or bimodal, with the maxima at k and n — k for some < A; < n/2, but not 
trimodal etc. We will assume this to be true in this paper but a formal proof is still lacking. 
Note that if we allow negative values of p the sequence can have a local maximum at every 
alternate index k. 

We write /(n) ~ g{n) to denote that f{n)/g{n) — )■ 1 as n — )■ oo. Analogously f{n) oc g{n) 
denotes that f{n)/g{n) — > A, for some non-zero real number A, as n — t- oo. 
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III. THE p,g-BINOMIAL DISTRIBUTION 



Now we are ready to introduce the notation 



fe=0 



(27) 



P>9 



and define the p, g-binomial probabihty function 



Pp,,(n, k) 



(28) 



The reader should here observe that the sum of the coefficients has, to the best of our 
knowledge, no simpler expression in the general case. Neither do the sum of the g-binomial 
coefficients have a simpler expression that we are aware of. Compare this with the case of 
the standard binomial coefficients for which the sum is simply 2". 

Having made our assumption of unimodality/bimodality we can now set up a simple 
computational scheme to find values of p given q. First we need to define a highly useful 
quantity; the ratio between two coefficients 





n 






_k 




i_ 






n 








_k_ 


p,g 



0<i<k<n/2 



(29) 



In the special case when £ — 1 we are looking at two consecutive coefficients. The ratio then 
becomes 



Rp,q {n, k, 1) = 





n 










1 


p,<i 




n 








A 


p,g 



pk _ qk 



pU—k+l qn—k+1 



(30) 



P 



~{n-2k+l) 



{q/pT 



1 - (q/p) 



n-fe+l 



-(n-2fe+l) 



1 - ip/<ir 



1 - {p/qr-'+' 

Lemma III.l. Let k,n be positive integers such that 1 < k < n — k. If x > 1 then 



< 



1 



X 



< 



n-k+i n-k + 1 



8 



and ifO<x<l then 



Proof. Obviously we have 



k 1 — 

< r-r < 1 



0<-I^<l 



]^ rjr.n—k-\-\ 

for all a; > so we proceed to the other inequalities instead. Define 

^ 1 + X + • • • + 

so that 

1 1 + a; + h x^~^ \ — x^ 



If X > 1 then 



rfM j^^k ^ Y X 



1 + y l + x-\ h 1 - x^-'^+i 

{n-2k+l)x'' n-2k + l n-2k + l 



1 -)- 1 _|_ . . . _|_ 2;. 



and then 



kx 



1 A; 

< 



k 



■ X > 



k 



1 + y n — k + 1 
If < X < 1, then completely analogously 

n-2k+ 1 n-2k+ 1 
y S ; X < 



and thus 



and the lemma follows. 



k k 
k 



> 



1 + y n — k + 1 



□ 



Theorem III. 2. Let k and n he positive integers such that k < n — k. For p > q > we 

have 



< 





n 










1_ 


p,q 




n 








>. 


p,g 



^ -(n-2fe+l) 



k 



n-k+1 



while for < p < q we have 



-(n-2fe+l) 



k 



n-k + 1 



< 





n 






> 




1 






n 










P>9 
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Given q and a coefficient ratio r at coefficient k we can now find the correct p through a 
simple iteration scheme: If 

< ^ < 5-(-2^'+l) ^ 



n — k + 1 

then p > g > and we use 



p ^ (g'^ - r g"-'+i + r p^-'^+i) ^ (31) 



which is obtained from setting f l30|) to r. As start value of p we may use a number slightly 
larger than q. If 



^-V"-^fv-t-±; 

n — 

then < p < g and we use 



^-yn-..^., < r < g- 

n — k + 1 



(Y)^ \ n — fc+l 

P_ + qn-k+i _<i_\ (32) 

and use as the starting value for p. To prove that these iteration schemes actually converge 
one would have to show that their derivatives with respect to p is at most 1 using the start 
value. Since we have no such proof we will leave it at that and just claim that they are 
practical. 

Given a p, g-binomial distribution, or a distribution that we wish to approximate by 
a p, g-binomial distribution, can we find the pair p, g > that generated it such that the 
distribution of p, g-binomial coefficients have the correct probability and ratio r at coefficient 
kl An iteration method again solves this problem practically under the assumption p > q. 
Suppose this input distribution has the probabilitites P(0),P(1) ...,F{n). Let k, 
r = F {k — 1) /F (k) and an e be given as input parameters. 

Algorithm p, g-Find 

1 

1. Assign gmin <- and gmax ^ ( r-(n-fc+i) ) " 

2. g ^ (gmin + gmax)/2 

3. Compute the corresponding p as in the method above. 

4. If Fp^g (n, k) <F {k) then gmin ^ g, otherwise gmax ^ I- 

5. If gmax — fein < c then exit loop, otherwise jump to step 2. 
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This method seems to work best when P (k) is one of the maximum probabihties. However, 
if the distribution is unimodal so that the maximum probabihty is at k = n/2, then the 
scheme will depend heavily on the quality of r. On the other hand, if the distribution is 
bimodal then this problem goes away and we may simply set r = 1, unless n is too small. 

It is implied, though we do not have a proof, that increasing q while keeping k and r 
fixed also increases the probability Fp^g (n, k). It actually increases until p = q which then 
constitutes an interesting limit case, which we will deal with in section IIVI 



IV. THE SPECIAL CASE p 



We will extend the definition of the p, g-binomial coefficients in (HM to include also the 
limiting case when p = g as in f[T5]) . Thus we will define 



k (n— fc) 



<?><? 



n 



(33) 



Let us look particularly at the point q where the coefficient ratio is 1 at k = n/2, i.e. 
Rq,q{n, n/2, 1) = 1. Also, henceforth we will assume that n is even to simplify some 
calculations. 

Lemma IV. 1. For q = we have Rq^q (n, n/2, 1) = 1. 
Proof. In the case when p = q the ratio is 



Rq,q{n, k, l)=g-(«-2'=+l) 



and for k = n/2 we have 



Rq,q {n, n/2, 1) = g ^ 



n — k + 1 



n 



n + 2 



Setting this to 1 gives the lemma. 



□ 



What is the sum of the coefficients at this point? To answer this we compare the middle 
coefficient with a coefficient situated at some carefully chosen distance from the middle. 
How big is the middle coefficient? Note first that 



n 
n/2^ 



q,q 



n 
n/2 



(34) 



11 



Lemma IV. 2. For q = we have 

^ 71+2 



n 



q,q 



2e / 2 



meaning that the quotient between the left- and right-hand side goes to 1 as n — )■ cxd. 
The proof follows from an easy application of the identity 



x\ " ^ / ^ 
1 + -) =e^ l- — + — r + — 



(35) 



n/ \ 2n 3 77.2 

and we leave it to the reader. A somewhat more involved application of (135|) is the following 
lemma 

Lemma IV. 3. Let x be some real number. For q = we have 

^ 71+2 

Rq,q {n, n/2, xn^^'^) ~ exp (^-^ . 
This allows us to give the exact order of the sum. 
Theorem IV. 4. For q = we have 

^ 71+2 



■ X' 



r(i/4) 31/4^1/4 



TT 



Proof. The calculations goes as follows though we leave out some details. 



n 
n/2_ 



n/2 

E 

9.9 k=-n/2 





n 






n 


/2 + k_ 


q,q 




n 








n/2_ 


9.9 



n 



3/4 



n 



-oo 



^3/4 /2e / 2 y 3V^r(v4) 
V nn VVe/ 2^2 
r(i/4) 3V"nV^ 

where the factor n^^^ in front of the integral comes from the change of variables k = xr?l'^. 

□ 

This is the sum at just a single point where the distribution becomes flat in the middle 
region. We can do even better if we allow ourselves to move around in the vicinity of this 
point. 
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Lemma IV. 5. Let 



for some real number a. Then 



n 



n a 



n + 2 n^/^ 



<?,<? 



2e / 2 



1 ^ I exp 



a \ n 



This can be verified using (!35|) as can the following lemma. 
Lemma IV. 6. Let a and x be real numbers. For q 

,3/4 



-1-^ H — TH we have 



Rg,g {n, n/2, xn^''^) ~ exp ^- 



4 

2 _ _ 4 

3 



We now have the resources to estimate the sum of the coefficients for a whole spectrum 
of values of q near n/(n + 2). The next theorem can be shown using the same technique as 
Theorem IIV.4t though the result gets slightly more complicated due to the integral on the 
right hand side in the previous lemma. 

Theorem IV. 7. Let q = H — For a > the asymptotic order of'^q^q (n) is 



n 



exp 



a \/n 3 



4 V TT \y/e J '""^ V 4 
For a < the asymptotic order of (n) is 



+ 



32 



^1/4 (^) 



n 



1/4 ^ 

, a/— 3 a e TT I I exp 
4 Vv^,' 



a ^ 3 



(/l/4(#)+/-l/4(#)) 



4 32 

Here Jq, (x) and -K^q (i^) denote the modified Bessel functions of the first and second kind 
respectively. The mth moment is simpler to express using an integral formulation. 



n/2 

E 1*1 

k=-n/2 

+00 



n 



q,q 



3 m+3 

n 4 



n 




' 1 im / 


n/2 




X , exp 1 







-ax 



^ - x'' 1 dx 



(36) 
(37) 



and the asymptotic behaviour of the middle coefficient is given by lemma IIV.5I 

The same technique allows us to repeat this for points farther away from the critical point 
q = n/{n + 2). If we increase q by a/n then the coefficients get sharply concentrated in the 
middle like that of standard binomial coefficients. Again ([35]) to the rescue. 
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Lemma IV. 8. Let q — ^ where a> 0. Then 



n+2 

n 



2e / 2 \" /an 
.^/2J V ^ W^J ''''^ V~ ~ ^ 

Note that above, when the coefficients had a rather wide distribution, we examined their 
behaviour at xn^^^ from the middle. Under our current assumption of q the distribution 
gets more sharply concentrated around the middle (basically they become gaussian), thus 
we study their behaviour at x y/n from the middle. 

Lemma IV.9. Let q ^ ^ and a > 0. Then 

Rq,q (n, n/2, X y/n) ~ exp (— ax^) 
Theorem IV. 10. Let q = ^ where a > 0. Then 



n 




_n/2 


/ 



+ 00 

y exp (— ax^) dx 




le I \ an a 



a \y/e ) \ 4 

If we decrease a below zero the sequence becomes sharply bimodal, with all its mass 
concentrated around two peaks. We can of course connect the position of the peaks with 
the parameter a. Suppose that we want one of the peaks to have its maximum located at k 
and A; — 1, that is, let the ratio here be 1. 

Lemma IV. 11. Given a number /i such < < 1 let 

atanh /j, 



a=2 [1- 
then 

lim Rg g {n, A;, 1) = 1 

n—>-co 

Proof. A simple calculation shows that the limit of the ratio is 

lim Rg g {n, k, 1) = lltif: exp (/x (a - 2)) 
Setting the limit to 1 and solving the equation gives the lemma. □ 
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We continue as before and estimate the growth rate of the peak coefficient. The result 
(and the proof) is somewhat more comphcated but follows from an application of f l5^ . 



Lemma IV. 12. Let fi and a be defined as in lemma UV. 1 1\ and set q = Th 



n 



Li (i + /i)J 



V2 



cxp 



atanh /i j + ^ (2 /i— atanh fi) atanh /i 



en 



A/vrn (1 - /i2) 



Again we take x y/n steps away from the peak and find the shape of the distribution. 



Lemma IV. 13. Let fi and a be defined as in lemma UV. 1 1\ and set q = Th 



en 



Rq,q \ (1 + /^) , a; ^/n^ ~ exp |2 a;^ {jH^^ ^ 



atanh /i 



Finally we get the sum by multiplying the integral with the peak coefficient and 2 y/n, 
where the factor 2 is due to that we have two peaks. 



Theorem IV. 14. Let fi and a be defined as in lemma UV. 1 1\ and set q = Th 



n 



+ 00 



'1,Q 



exp <2x 



1 atanh /i\ , , 
; + 1 dx 



en 



V2 



TT n 



n 

I (1 + /^). 



q,q 



yU + (/i^ — 1) atanh /i 



where the growth rate of the peak coefficient is that of lemma \IV.12i 



V. POCHHAMMER BOUNDS 

In this section we set up bounds for the g-Pochhammer function and use them for giving 
bounds of quotients between g-binomial coefficients. Basically we mimic the upper bound 
in^ andifi but we extend it to obtain a lower bound as well. They are very useful bounds so 
we will do this in some detail though everything is based on standard elementary methods. 
First we need the integral estimate of a sum. Let f{x) be a continuous, positive, decreasing 
function on the interval m < x < n + 1 where m and n are integers. Then 

n+l ^ n 

J fix) dx<J2 /(^) ^ /(^) + / /(^) (38) 

_ k=m „ 
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Recall that the dilogarithm is defined as 

oo 

Li2 (x) = 

n=l 

Now let < a < 1 and < g < 1 and note that 



log(l-t) 



dt 



(39) 



n-l 



- log (a; g)^ = ^ - log (l - a q'') 



(40) 



k=0 



k=l 



Note also that — log(l — aq^) is a positive and decreasing function for a; > 0. Take the 
series expansion 

oo 

-log(l-a;) = ^y 

k=l 

SO that 

-log(l-ag^) = E^"'^ 

Integration gives 



k 



(41) 



(42) 



A:=l 

g (7 I ^ A;2 ^ 



log q 



\ogq 



,fc=l k=l 

Together with the integral estimates above we have 

Li2 (ag") - Li2 (a) 



k=l 

Li2 (ag") - Li2 (a) 

logg 



(43) 
(44) 



logg 



< - log (a; g)„ < - log(l - a) + 



Li2 (ag"-i) - Li2 (a) 



logg 



(45) 



Reversing the signs and taking exponentials we finally obtain 

'Li2 (a) - Li2 (aq" 



(a; g)„ > (1 - a) exp 
(a; g)„ < exp 



logg 

Li2 (a) - Li2 {ag" 



logg 



(46) 
(47) 



Now we turn to the g-binomial coefficients. Let < i < k < n/2 and use (|5]) to note 
that the ratio between coefficient k — i and coefficient k is 



Rg (n, k, 





n 






_k 




i_ 


g _ (g^ ^+1; g)^ 




n 




(g"-^+i; g)^ 




_k_ 







(48) 
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Using the bounds for the g-Pochhammer function we can now bound the ratio. For the 
upper bound of the ratio we take the quotient of the upper bound and the lower bound. 
The ratio Rq (n, /c, tj then has the upper bound 



exp 



' Li2(g''-^+^)+Li2(g"-'°+^)-Li2(g''+^)-Li2(t,"-fe+i) ' 

R, in, k, i) < ^ 1 _ (49) 



and, quite analogously, it has the lower bound 

H, (n, k. t)>{l- exp ( "'''""'^"'''-';':;]-"'''')-"''''"") ) (50) 

Using (IT7|) we can now obtain bounds for the quotients of p, g-binomial coefficients. Note 
simply that 

Rp,q {n, k, i) = ^^^^^_2k+i) v^l) 
and use the bounds from fHOll and fl50l). 



VI. CONTROLLING THE p, g-BINOMIAL DISTRIBUTION 

Choosing a for a given n such that Rp^q (n, /c, 1) = 1 defines a set of pairs p, g; an 
isocurve. If we choose a value of g, then what value should p have to result in a distribution 
of coefficients which has a peak at k (and n — k), i.e. with Rp q (n, k, 1) = 1? The iterative 
method flHTl) above produces the correct p for any given g but reveals no information on p. 
To obtain this we need to parameterize g properly and one way to do this is to set g = H-2;/n 
for some z < 0. This parameterization was also used in^ andiS- for computing giving upper 
bounds on g-binomial coefficients. 



A. Wide and flat distributions 



Let us begin with the particular distribution which has its peak in the middle, i.e. k = 
n/2. Due to symmetry we are interested only in the case p > q. Recall that in the special 
case p = q the distribution has its peak at n/2 when 

n , 2 4 8 

+ + (52) 



n + 2 n 
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which would correspond to z = —2. Let us also work under the assumption that p has the 
expansion 

p = i + yi + yi + ... (53) 

n 



+ ■■■ (57) 



and determine what values of yi,y2, ... we should have. First, according to fl30|l . we should 
have 

r,n/2 _ n/2 

Rv,, in, n/2, 1) = ^ J,^, _ = 1 (54) 

which we rewrite as 

pn/2+l _ pa/2 ^ ^n/2+1 _ ^n/2 (55) 

Setting p = l + yi/n + y2/n'^ + ■ ■ ■ and q= 1 + z/n and performing a series expansion using 
f l5^ . we find the right hand side to be 

z z^{?,z + 16) e^/^ 

~ 4^ ^ 967^3 ^ ■ ■ ■ ' 

and the first two terms of the left hand side are 

yiey^'^ {yl-Ay2-2yiy2) ey^'^ 
n An^ 

We solve this term by term. First 

y^e^'/^ = ze'^^ (58) 

has the solution yi = 2w where 

w = w(^^ j (59) 

Here s = W (x) is the Lambert function solving se^ = x. Note that for z > —2 we have 
yi = z but this is not the case for z < —2 where we have yi > z. The second coefficient is 

2(u,+ l) 

with w as before. We could go on and solve for y3,y4, . . . but, for the case in hand we 
actually only need yi. We will henceforth drop the subscript and refer to it as simply y. 

What shape does the distribution have at this particular p and g? What we are seeking 
is an expression for Rp^g {n, n/2, i) and this is where we start using (149|) . (l50i) and (15T1) . We 
define i = xr?!'^, just as we did in the case oi p = q. First we need an expression for q/p 
and we let 

p l + ^ + --- n ^ ' 
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where w is defined as in We are now ready to compute the hmits of the upper bound 
din]) and lower bound f lSUl) . It turns out that the limits of these bounds coincide and we 
receive 

(n, n/2, X v?"') ~ exp f ^—^ ^— \ (62) 

Note that the special case z = —2 corresponds to p = q and gives the coefficient —4/3 of x^. 
The calculations were performed with Mathematica and are much to unwieldy to fit in this 
paper. We have prepared a Mathematica notebook that performs the calculations step by 
step, using some practical transformation rules. What we compute is actually the limit of the 
logarithm of the upper and lower bound. The steps are as follows; compute series expansions 
of the different powers of r = g/p, use them inside the dilogarithms and then compute their 
series expansions, add the dilogarithms and the series expansions of the logarithms of the 
remaining factors. Some transformations of this expression helps Mathematica to take the 
limit that gives the result. 

Figure [T] demonstrates how the asymptotic ratio is achieved with increasing n. It shows 
Rp,q {n, n/2, xn^/^) ior z = —9 at y = 2w = —0.10539 . . . and the asymptotic ratio is given 
by (ES]), that is, e-^-^^^^^^^ . The red curve is the asymptote and the blue curves are for 
finite n where the curves for larger n are closer to the asymptote. 



B. Wide and double-peaked distributions 

Changing p only slightly, say on the order of 1/n^/^, allows us to move around in the 
region where the distribution is wide. The end result is that with q = 1 + z/n and p = 

I + 2w/n + a/ r?!'^ we get 



„ / , Q/Ax (wz(2w\z^x^ az{WvS)x^\ 

Rp,, n, n/2, xn'/^) ~ exp ^— + ^— 

^ ^ \ 6 2w — z J 



(63) 



The ratio expression gets nicer if we change the coefficient of 1/n^/^. We suggest the following 
parametrization instead; with g = 1 + z/n let 

2w aw iz"^ — Aw"^) 

P=^ + — + . 3/2 64 

n S{1 + wj n-^'^ 



and we receive 



Rp,, {n, n/2, x n'/') ~ exp (^!!ili^^±A [x'~2ax')^ (65) 
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FIG. 1: Rp^q [n, n/2, xv?f^^ (blue curves) versus x at z = —9 and y = 2w = —0.105391 for n = 2^, 
2^, 2^° and 2^"^. The red curve is the asymptote e'^-'^^^'^^^''^ . 

This puts the maximum a.t x = ±y/a for a > and at x = for a < 0. 

Figure [2] works hke figure [1] but here with the parameter a set to 1. It shows 
Rp,q {n, n/2, xn^/^) for 2; = —9 with p = I — 0.105391/n— 1.50171/n'^/^ and the asymptotic 
ratio is given by dnS]), that is, e-o-7i97i8 (x4-2x2)_ r^y^^ 

curve is the asymptote and the blue 
curves are for finite n where the larger n are closer to the asymptote. Had we set a < the 
distributions would still be wide but with a single peak in the middle. 

C. Peakish distributions 



If we instead choose a middle ratio of 1+a/n then we receive a sharply peaked distribution. 
Say that we want 

p"/2 _ qn/2 ^ 

Rp^q (n, n/2, 1) = — ^- -j^ = 1 + - (66) 
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FIG. 2: Rp^q (n, n/2, xn^/^) (blue curves) versus x at z = —9 and a = 1 so that q = 1 — 9/n and 
p = 1 -0.105391/n- 1.50171/^3/2 for n = 2^, 2*^, 2^0, 2^2 and 21"^. The red curve is the asymptote 

g-0.719718(a;*-2a;2) 

for p = 1 + Hi/n + ■ ■ ■ and q = 1 + z/n a.s before. Note that a = —2 corresponds to a 
binomial distribution, so we are usually interested in the case —2 < a < 0. After expanding 
the equation we receive as before, from the first term of the left and right hand side, the 
equation 

(a + yi) = (^a + z) e'/^ (67) 
which has the solution yi = 2w — a where 



2 

The second equation, which is slightly longer, has the solution 

_ {a-2w + z) {a'^ {w + 2) - 2aw'^ -w z{2w + z)) 

- 2(w + l)(a + z) ^^^^ 
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but we do not really need it at this moment. This gives a distribution with a width of the 
order ^Jn. Computing the limit ratio gives us 

(n, n/2, X v^) ~ e"^' (70) 

thus giving us an essentially gaussian distribution. Note that it does not depend on z in 
its current form. Of course, we expect other terms to depend on z but these vanish when 
n — oo. 



D. Two separate peaks 

Suppose that we want the peaks located outside the middle. Defining k = | (1 + /i) for 
< < 1 means that we move the peaks out from the middle and that we move to another 
isocurve. We keep g = 1 + z/ra and p = 1 + yxjn + y^jTi^ + ■ ■ ■ and solve Ryq (n, A;, 1) = 1, 
i.e. 

With /c = I (1 + yu) we use (l35ll on both sides and find the equation 

exp (I (1 - /i)) - exp (I (1 + /.)) = exp (1 - - exp (| (1 + (72) 
which we rewrite as 

eS'i/2sinh^ = e^/2gi^j^/i£ (73) 

At this point it would be appropriate to define the function s = Qfj_ [x), for < < 1, as 
the maximum solution to the equation 

X = sinh (s /i) (74) 

Note here that the function sinh(yux) has a minimum at (— atanh/x)//i for < /i < 1 
and, due to symmetry, a maximum at the same point for — 1 < yU < 0. The function fl 
returns a value in the interval 

atanh/i 

< s < (75) 

H 

The solution sought in our equation is thus yi = 2w where 



w 



% sinh ^) (76) 
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We should mention that (x) is a natural extension of W (x). In fact, if we let w = W (x e^) 
then 

fi^(e^ sinh/ix) =w + — ^— + (77) 

b[l + w) 

giving a good approximation for small values of fi. The reader may recall (1601) above for 
comparison. 

Having computed p = 1 + yi/n we compute the upper and lower bounds of the ratio as 
before. Unfortunately, the ratio has the rather ghastly expression 



Rp,q ^ (1 + xVn) 



r x2 / e^''"+2 {e^^ - e') {z - 2w) \ 1 

where w is defined by ( 176|) . 

In figure [3] we show how the finite cases approach their asymptote for /i = 1 /3 and 
z = —9. This gives p = 1 — 0.153208/n and the asymptotic ratio is given by (178|) . that is, 
g-o.io355ix2^ The red curve is the asymptote and the blue curves are for finite n where the 
larger n are closer to the asymptote. 



VII. MOVING ALONG A DIAGONAL 

With q = 1 + z/n and using the expressions given by flM|) and (jHSD , we can move around 
in the region where Rp^g {n, n/2, 1) is very close to 1. It is understood here that a change 
in the parameter a only changes p, i.e. we only move in the p-direction and keep q fixed. 
Let us instead say that we want to move, ever so slightly, in both the p- and g-direction. If 
the movement is on the order we can translate it to a different parameter a in the 

expression given by flM|) . 

Given a z in q = 1 + z/n the starting point for p is p = 1+y/n with y = 2w and w given 
by ( 159|) . Let us think of the starting coordinate as {yo, zq). We wish to move from this point 
sj ^fn steps in the direction t, where the limit t — ?■ oo corresponds to the case where we only 
move in the direction y (or p) . The new point is now 

and note here that \/ {yi — y^Y + {zi — zqY = |s| / ^/n. 
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FIG. 3: Rp^q (n, n/2 (1 + x ^/n) (blue curves) versus j; at /u = 1/3 and z = —9 so that q = 1 — 9/n 
and p = 1 - 0.153208/n for n = 2^, 2*^, 2^°, 2^^ 2^^ and 2^^. The red curve is the asymptote 



-0.103551 



Having moved in the ^-direction from zo to zi we now compute the point ?/2 so that 
(?/2, Zi) stays on the isocurve with Rp g {n, n/2, 1) = 1. A small change in z requires only a 
small change in y to preserve this property. 

With wq = W (^zq e^"/^/2) we have the series expansion 



^zq + x) exp{{zo + x)/2)\ xwo{zq + 2) 

W I z ] = Wo + 7 + 



2 y " ' 2zoK + l) 

giving an approximation for small x. We want to set 



X 



and thus we have 



1/2 = 1/0 + 



swo (zo + 2) 



zo (wo + 1) VI + 



n 



(80) 



(81) 



(82) 
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The difference yi — 1/2 is what we want: 



st swo{zo + 2) 
Vi — V2 = , , = (83) 



ym^y^ V ^0(^^0 + 1) 



3 [wq + 1) x/n 



(86) 



Now solve 



where the right hand side is extracted from flMl) . This gives 

3 s 1 2:0 (wo + 1) - Wo (-20 + 2) 



which then gives the ratio in (jHS])- 

As an example, we take the case p = q which corresponds to t = 1 and zq —2~. Taking 
the limit of the expression for a gives 

3s 

lim a = 1= (87) 



which corresponds to the coefficient ratio 



Rp,q {n, n/2, xn^^"^) ~ exp i^'^ ~ 2ax^)^ = exp ^- 



-x^-^xM (88) 
3 y/2 ' ^ ^ 



Compare this with lemma [TV. 6 j the parameter a in the lemma corresponds to taking s = a \/2 
steps in direction t = 1 and this gives the same coefficient ratio as in the lemma. 



VIII. MOMENTS 

Once we have the ratios (|^ . fl^ . flTUl) . (1751) it is an easy task to compute moments of 
the distributions. Let us do this for the most interesting case of (165|) . First, to make the 
notation somewhat simpler, denote = $(z) = —wz {2w + z)/6, i.e. (p > 0, with w as in 
f l59|l . so that 

Rp^g {n, n/2, xn^/^) ~ g-'^^''-^'^"') (89) 
where a is defined by (^^. We use the notation 

) (90) 
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for the mth moment of the probabihty distribution of k = 0, . . . ,n. Define also 

+00 



J Ixl'" exp (— (x^ — 2 ax^)) dx, m>0 



so that the mth moment becomes = Qm/Qo- For m = we have 



exp (0a^/2) (/1/4 ((/-aVa) + (^aVa)) , for a > 



= < 



^ exp(0aV2) iri/4(<^aV2) 



for a < 
for a = 



so that 



Next, for m = 1 



n/2 



k=-n/2 



270 



exp (0 a^) ^1 + erf (^a a/0 



where — 1 < erf fx) < 1 is the error function, and thus 



0-1 



n 



3/4 Ql 
Qo 



In general we have for m > that 



n 



k-- 
2 



n 



3m/4 

^0 



where is given by 



2 (r(in±i) iFi(in±i,i,<^„2) + 2av^r(i2±3) iFi(i2±2,|,0a2)) , a>0 
<i^SS7?r(-±i) U(^,i,0a^), a<0 

a = 



2 V 4 7 ) 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 



(97) 



Here iFi(a,6,c) and U(a,fe,c) denote the confluent hypergeometric functions of the first and 
second kind respectively. If we want to compute cumulant ratios we first need moment ratios 
which of course is easy now. For example, in the case of a = we have 



02_ QqQ2 



72 = 1.4142... 



?A - . ^0 g4 _ r(i/4) 



2.1884. . . 



(JQ qI Qq 3r(l/4) 

~ — 5— = — = O.OOOO - 



(98) 
(99) 

(100) 



These moment ratios are the same as obtained in the 5-dimensional Ising model, see e.g 



11 
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IX. FINE-TUNING THE EXPONENTS 



Note (1^ and, as before, keep $(2;) = —wz {2w + z)lQ where w is defined by (15^ . Recall 
that the first and second absolute moments obtained from (1971) for a = are 



00 r(i/4) $(^)V4 $.^xi/4 ^ ^ 



■?/9 f'o 0/0 TT v2 1 ri'^/^ 
^2 ~ n^/^ ^ ^ ^3/2 V ^ 

f?o r(i/4)' y^^z) 



The argument z is allowed to depend on n but probably not to a high order. At this point 
it is not clear how z may depend on n for the calculations leading to fl65p and ^{z) to work. 
We will assume, for the moment (see the end of this section), that the expressions for the 
moments of fl97|) are valid when z = O (logra). 
The series expansion of W (x) is 

w = W (x) = X - + — — + ■■■ (103) 

With X = e^l"^ z/2, and note that z is negative, we have 

,,,, , ^^Jze'l'^\ ^2 z 3^3 3^/2 

^ = W(x)=w(^^J =_ — + + ■ ■ • (104) 

so that 

, — 2; w (2 ti? + 2;) 2;^ e^/^ z^e^ z^ , . 

$(z) = ^ = \ \ (105) 

^ ' 6 12 12 24 ^ ' 

Set z = Xq + Xi logn + A2 log log n + A3 log log log n with Ai, A2 < and focus on the first 
term of ffTU3]) . 

-1 3 
$(-2) ~ — (Ao + Ai logn + A2 log log n + A3 log log logn) 

We are interested in two special cases. First choose Ai < 0, A2 = —6 and A3 = 0. This gives 
Combining this with f llOip we receive 

27 



and 

' r(i/4)^ (-A,)'/2e^./' 
Had we let A2 = 0, instead of A2 = —6, then we would have ended up with a factor log^^^ n in 
the denominator of fllOSp and a factor log^^'^ra in the denominator of f llOOp . For the second 
case we choose Ai = 0, A2 < and A3 = —6. We get 



This together with f llOip gives us 



31/^ 727 ^3/4 log- 



n 



CTi ~ (111) 

' r(i/4) (-A2)3/4e^o/8 

and 

These expressions obviously converge extremely slowly and are probably not of any use for 
n that might occur in practical situations. 

We have managed to verify fllOSp and (11091) by using the method described in section \V1\ 
again using Mathematica, only in the special case z = — \ogn, i.e. Ai = —1, Aq = A2 = A3 = 
0. We could then confirm that 

2; ^7/8 \ ( —X^ 



) ^ \J2 ) ^^^^^ 

Computing the moments of this distribution produces the same result as setting Ai = — 1 
and Ao = A2 = A3 = in fll06p and then computing the moments in the same the way we 
obtained fllOSp and fllOOp . A more general computation seems not to be within reach with 
our current set of tools though. To conclude this section we note that the exponent of n in 
(11 13p is 3/4 — Ai/8. For this exponent to stay less than one we thus need Ai > —2, giving 
us a bound on z. 



X. THE ISING MODEL 

A state r on a graph G is a function from the set of vertices to {±1}. There are thus 2"' 
states for a graph on n vertices. We define the energy of a state as E{t) = '^■j TiTj where 
the sum is taken over all edges ij of G. The magnetisation is defined as M(r) = J2i '^i with 

28 



the sum taken over all vertices i of G. Note that —n < M < n and it only takes every 
alternate value, i.e. M e {— n, —n + 2, —n + 4, . . . , n — 4, n — 2, n}. We will often need to 
refer to it in terms of how many negative spins the state has. If k spins are negative then 
M = n-2k. 

The partition function of the Ising model is defined for any graph G as 

r E,M 

The coefficients a(i?, M) then are defined as the number of states with energy E and mag- 
netisation M. Denote the number of states at energy E by a{E) — ci{E, M). Note that 
the number of states at magnetisation M is just (^), where k — {n — M)/2 is the number 
of negative spins. Let also denote the terms of Z with magnetisation M = n — 2 A; for a 
graph on n vertices, so that Z = Zq + Z2 + • • • + Z„, i.e. Z^ are the terms corresponding to 
k negative spins. 

If we evaluate the partition function in x = and y = with K the dimension- 
less coupling, or inverse temperature J/ksT, and H — h/ksT as the dimensionless exter- 
nal magnetic field, we obtain the physical partition function denoted Z — Z{G; K, H) — 
Z(G; e^, e^), though we are usually interested only in the case when if = (or y = 1). 
Analogously, we write Z — Zq + • ■ ■ + Zn- The dimensionless and normalised free energy 
is defined as = (logZ) jn. Prom the derivatives of the free energy we can now obtain 
other physical quantities such as the internal energy dT jdK and the specific heat d^J^/dK'^ 
though we shall not be needing the latter for this investigation. 

We assume the Boltzmann distribution on the states so that (with H — {)) the probability 
for state r is 

KE{t) 

P(r) = ^ (115) 
We have then especially that the probability for energy E is 

HE)^^-^ (116) 
and the probability for magnetisation M is 



E 

where M — n — 2k. 
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Denote by K* the coupling where Zn/2-1 = Zn/2 = Zn/2+i- This couphng will corre- 
spond to the relation Rp^g{n, n/2, 1) = 1 for some choice of p,q. We define the (sponta- 
neous) normalised magnetisation /2 = {\M\) /n and the (spontaneous) susceptibility x = 
Var(|M|)/n = ((M2) - (|M|)^) /n. The pure susceptibility is simply x = Var (M) /n = 
(M^) = Aa2/n. Since M = n — 2k we thus have /2 = 2ai/n and x = 4 ((T2 — cyD/n. Recall 
the traditional finite-size scaling laws which claim that in the critical region, i.e. near Kc, 
Ji oc L~^/'^ and x ^ L'^^'^ ■ Being near Kf. means that \K — oc L~^/^ and we especially 
expect K* to belong to this region. Though the high- and low-temperature exponents may 
or may not be equal for three dimensions, see^ for an in-depth numerical investigation of 
this matter, the details of these exponents are not important for our present investigation. 
What matters is that there are exponents that guide the growth of e.g. the susceptibility 
near K^. 



A. The complete graph 

For a complete graph, denoted Kn on n vertices and (2) edges the partition function is 
easy to compute. Suppose k of the vertices are assigned spin —1 and the other n — k have 
spin +1. The magnetisation is obviously M = n — 2 k and the energy is 

E= + (''~^] -kin-k)= ("f] - 2kin-k) (118) 



2J \ 2 J ' V2 
The partition function is then 

k (n—k) / 1 \ 



fc=0 

and with y = 1 we have 



n 



ZiK^; X, 1) = xi"^) J2 



k=0 



= xUJ ^^J^ ,^ (n) (120) 

1,1 



where q = Thus we have 



n 



Z{Kn; K, 0) = exp | (n) (121) 

where q = exp{—2K). Obviously we have 

F{M = n-2k)= Fq^q (n, k) (122) 
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Since we have defined the critical temperature as the K = K* where the middle ratio is 
1, i.e. P (M = -2) = P (M = 0) = P (M = +2) then this corresponds to the point where 
Rq,q {n, n/2, 1) = 1, which takes place at q = n/{n + 2) as we saw in lemma llV.li Thus 
K* = i log (1 + I) for Kn. 

In short, the partition function and the magnetisation distribution for Kn can be ex- 
pressed in terms of p, g-binomial coefficients. Does this hold for all graphs? No. In fact, it 
seems to only be true for Kn- However, it does seem to hold asymptotically as the order of 
the graphs increase, for some interesting classes of graphs. The precise formulation of such 
a statement remains and falls outside this paper. 



B. The average graph 



Let us compute the sum of all partition functions taken over all graphs on n vertices. 



= Z(G; x, y) 

GCK„ 



i=0 



("D 



E 

i=0 



n 



E 

k=0 



E 

^=0 



X 



i {n-—i)\ 

1 \ I 2 
1 + - 

X 



n 

{1 + xp y^Y. 



i=0 



n 



1 \' I 



X 



i (n—i) 



(123) 

(124) 

(125) 
(126) 



and for y = 1 we have 



Z„(x,l) = (1 + x)^^^ Y 



n / \ / -i \ i (n—i) 



i=0 
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(127) 
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where q = 1/x so that K* = log(l + 2/r;,). Again we have P (M = n — 2k) = P^^g (n, k). The 
mean magnetisation distribution can then be modelled by p, g-binomial coefficients, though 
with p = q, just as for the complete graph. 
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C. The complete bipartite graph 



What about K^^^, i.e. the complete bipartite graph on n = u + v vertices? Now the 
partition function is 

Z{K^y, X, y) = (129) 

^^^^^ yU+v-2i-2j ^ij+{u-i)iv-j)-iiv-j)-jiu-'i) ^ ^-j^gg^ 



X 

i=0 j=0 



which for y = 1 gives us 



Z(ir„,,; X, 1) = (132) 

EE :K Ui) - («3) 



i=0 j=0 



k=0 I 

which defines the partial sums for ?/ = 1 as 



EE ''''' 



Data suggests that 



¥{M = n- 2k) = ^ Pp,g (n, fc) (136) 



given an appropriate choice of p and q and for a rather wide range of temperatures. It 
does however not seem to hold if u differ from v. In the left panel of figure H] we show 
a sample of magnetisation distributions together with fitted p, g-distributions for a -^'32,32- 
To find the appropriate p and q we used the method described in the p, g-find algorithm 
in section lllli The fit is excellent. The right panel of figure H] shows y = n{j) — \) versus 
2; = n (g — 1) for a range of temperatures and for complete bipartite graphs of different sizes. 
High temperatures, K = ^ gives p = g = 1, i.e. ?/ = z = in the upper right corner. As 
the temperature decreases, i.e. with increasing we move along the curves. The points 
are where the distribution is exactly flat, i.e. the inverse temperature K* where the middle 
probabilities are equal. We have no exact closed form expression for K* but one can show 
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FIG. 4: Left: magnetisation distributions for ^32,32 (points) and fitted Fp^g (64, k) (lines) vs /c — 32 
for four different temperatures. Right: y = n{p — 1) versus z = n{q — 1) for Kn/2,n/2y with 
n = 24,32,48,64,96,128,192 (downwards). Points represent K* . Low temperatures in lower left 
corner. 

that the series expansion of this inverse temperature for Kn/2,n/2, i-e. a total of n vertices is 

K* ^ - 3827 

The calculations behind this are rather long and were done with Mathematica. Compare 
this with the expansion in the previous subsection for the complete graph on n vertices, Kn, 
which begins K* = 1/n — l/n"^ + 4/3n^ + ■ ■ ■ . In the lower left corner the distribution has 
P (M = n) = F (M = n — 2), i.e. at K = logn/n. The points in the right panel of figure H] 
should approach z = y = —2. One can show that the distribution of a balanced bipartite 
graph Kn/2^n/2 has the shape 

R {n, n/2, x n^'^) ~ exp (138) 

at K* , just like the complete graph on n vertices. Since $(—2) = 4/3 we then assume that 
z, and thus also y, will approach —2. 



D. The free energy 

If the magnetisations were indeed an exact p, g-binomial distribution then we could also 
express the free energy as 

HG-^K) = '^ + ^^^^^^ (139) 

n n 
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for a graph on n vertices and m edges and it is here imphed that p and q depend on K. 
Why this expression? Note that Zq = a{m, n) = e^"^ and thus we have 



Pp,9 (n, 0) 



1 



in) 



in) 



Z 



Z 



(140) 



from which the result follows. Compare with fll2ip where this relation holds exactly. Ac- 
tually we expect f ll39p to be a good approximation for K near 0, where the distributions 
are close to binomial, and for very high K where all the probability mass is concentrated 
on the extreme magnetisations. In the left plot of figure [5] we show the exactly computed 
free energy (red curve) for a complete bipartite graph on 16 + 16 vertices together with the 
p, g-approximation (11391) (points). The fit is indeed very good for the whole temperature 
range. Taking a derivative of the points with respect to K produces a good approximation 
to the internal energy as the right plot shows. 



dF/dK 

6r 



0.85 - 





FIG. 5: Left: free energy (red curve) compared to the formula (|139|) (points) for i^i6,i6- Right: 
internal energy (red curve) compared to the derivative of the points produced by ()139p (points) for 

-f^l6,16- 



XI. LATTICES 

The plots in figure H] are very representative for several graphs of interest. We intend to 
focus on the graphs that traditionally are studied in statistical physics; lattice graphs. We 
will take a look at the simple lattices in 1, 2, 3, 4 and 5 dimensions. Just to be clear, a 
ID-lattice is a cycle C„ on n vertices and it is 2-regular, i.e. 2 neighbours for each vertex. 
The 2-dimensional L x L-lattice is the cartesian product of two cycles on L vertices. The 
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product thus has n = L"^ vertices and it is 4-regular. The d- dimensional Lx Lx ■ ■ ■ x L-lattice 
is a product of d cycles on L vertices, thus having a total oi n = L'^ vertices. It is obviously 
2 (i-regular. Assuming finite-size scaling to hold then for a (i-dimensional lattice we have 

ai = nJji/2 (xnL-^'" = n {n^''^y^'" = n^-^l'^'' (141) 

and correspondingly for the second moment 

a2 = nx/4 oc nL^'" = n {n^''^y''' = n^+T/'^^ (142) 

Note also that for an r-regular triangle-free graphs we have P (M = n) = P (M = n — 2) 
when K=^-^. 

2 r 

A. ID-lattices 

For ID-lattices we can compute the coefficients a{E, M) exactly. It is an exercise to show 
that the number of states with k negative spins and i negative spin products (over the edges) 
is 

where M = n — 2 k and E = n — 2 i. The distribution of magnetisations do not behave in a 
way representative for lattices of higher dimension. However, for extremely low temperatures 
the probabilities P (M = —n) = P (M = n) will dominate the other probabilities. The two 
outermost probabilities, P(M = n) and P(M = n — 2), are equal when K = For 
the ID-lattice the distribution is here sharply unimodal, while for higher dimensions the 
distribution is bimodal and has its peaks at the extreme magnetisations. For K larger than 
(logn)/4 the distribution actually has three peaks, i.e. a local maximum at M = 0. For 
ID-lattices the p, g- approximation of the distribution thus breaks down beyond this K since 
it can not model a local maximum in the middle as well as peaks at the ends; they are at 
most bimodal. For K less than this point the p, g-distribution is a very good approximation. 
Figure [6] demonstrates this clearly; for the flattest distribution (low temperature) the fltted 
p, g-distribution starts to deviate from the actual distribution. In flgure [7] we plot y = 
n{p — 1) and z = n{q — 1) versus K for a range of different n. Clearly there is some limit 
curve here, though we have not established what the limit function is. In flgure M we see 
y = n{p — 1) versus z = n{q — 1) for different n. The right plot of flgure [S] shows the value 
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FIG. 6: Magnetisation distributions for C32 (points) and fitted Pp^g (n, k) (fines) vs A; — n/2 for 
several temperatures. 




FIG. 7: Left: y = n{p — 1) versus K for C„. Right: z = n{q — 1) versus K for C„. Both plots are 
for n = 16,32,64, 128 (larger cycles stretch farther to the right). 

at K that gives the maximum value of y. The fitted straight hne gives the hmit 0.1333, 
very close to 2/15. What about the values of y and z7 Indeed they converge beautifully 
as figure |9] indicates. The limit for y is about 1.010 and z approaches a value of —3.537. 
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FIG. 8: Left: y = n{p — 1) versus z = n{q — 1) for C„, with n = 8,16,32,64,128 with larger n 
extending farther to the left. Right: K giving the maximum y vs 1/n for Cn, n = 8, 12, 16, 24, 32, 
48, 64, 96, 128, 192, 256. 

Though we can not exactly solve what y and z should be at K = 2/15 we can at least see 
y z 
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256 64 32 24 16 12 

FIG. 9: Left: Maximum value of y = n{p — 1) versus 1/n for C„. Right: value of z (q — 1) versus 
1/n for Cn when y is at its maximum. In both cases n = 8, 12, 16, 24, 32, 48, 64, 96, 128, 192, 256. 

how y and z relate at this point. For an infinite 1-dimensional lattice we have that X = e^^, 
see e.g.—. The second moment then should behave as 



^2 ~ — = (144) 

Let i = k — n/2 and cr = y/d^. For high temperatures we expect £/a to be normally 
distributed and thus 

cr v2 TT 

The probability ratio is then 



2K 



R{n, n/2, i) 



¥(0) 



exp{-f/2a^) 



(146) 
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and for i = 1 this simplifies to 



2e 



-2K 



R{n, n/2, 1) = exp (-1/2 a^) = exp {-2e~^^/n) ~ 1 



n 



(147) 



Compare this with (166|) . We thus have a = —2 e~^^. Now ?/ and z are related as ?/ = 2 w — a 
where w is defined by fl68l) . If we set = 2/15 then a = —1.531857, and choosing z = —3.537 
indeed gives ns y = 1.01002. To actually solve 2; as a function of K seems harder though. 
However, numerical experimentation suggests that y and z for small K behave as 



where Ci ~ 6.164 and C2 ~ —10.33. 

Strangely, when the p, g-distribution fit the magnetisation distribution so well one might 
think that the free energy would be well approximated by f ll39p . This is not so. The 
p, g- approximation differs clearly from the asymptotic free energy, given by log (2 cosh if). 

B. 2D-lattices 

For the 2-dimensional lattices we can rely on exact data only for up to L = 16 
and they were computed according to the method We have sampled data for 

L = 32, 64, 128, 256, 512, collected with the methods described in^^ and^. These meth- 
ods gave us the energy distribution and then it is just a matter of combining this with the 
distribution of magnetisations for each given energy as described in^^. Figure [10] shows an 
example of some distributions for the 128 x 128-lattice together with their fitted p, g-binomial 
distributions. The fit is fairly good, but hardly excellent near K*. However, as the figure 
shows, at i^' = 0.4388 (i.e. for L = 128) the fit is practically spot on. For the lattices we 
have studied there is always one such temperature where the p, g-distribution fit particularly 
well. This point is located between K* and and is very close to, but not exactly equal 
to, the point where the susceptibility is at its maximum. 

Of course, for high temperatures (small K) and low temperatures (high K) the fit is typ- 
ically very good but in the high-temperature region the measured y and z are unfortunately 
extremely sensitive to noise. As we get closer to the critical region where the distribution 
becomes bimodal this problem goes away, even though the sampled distributions are more 



y{K) ^ ciVk + C2K 
z{K) ^ -ciVk + C2K 



(148) 



(149) 
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noisy there. Regarding the free energy it is well- fitted by ( I139p for low temperatures K > K* 
though less well for high temperatures K < K* . 
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FIG. 10: Magnetisation distributions for the 128 x 128-lattice (red) together with fitted p, (/-binomial 
distributions Pp,, (n, k) (blue) ys k - n/2 at K = 0.432, K* = 0.43467, K = 0.437, K = 0.4388 
and Kc = 0.44068 (downwards at y-axis). 



Recall from section HXl how the exponents of the moment growth rates could be computed 
if we allow z to depend on n. For the 2D-lattices it is known that /3 = l/8, 7 = 7/4 and 
z/ = 1, seei^,— anc^^. Thus the first moment ai should scale as n^^/^^ and the second moment 
(72 as n^^/®. From equation (11081) and (11091) this would be achieved by choosing Ai = —3/2, 
A2 = — 6 and A3 = 0. The left plot of Figure [TT] shows z versus logn at K* together with 
the curve 3 — 1.5 logn — 6 log logn. The constant Aq is chosen only to make the curve look 
plausibly near the points. The point for L = 512 deviate slightly but we suspect that noise 
in the sampled data explains this. With Aq = 3 the coefficient of n^^^^ obtained from (I109p 
would be 0.301 though the measured a2 divided by n^^^^ are closer 0.08. To get this we have 
to choose Aq ~ 8.3. In that case the convergence is extremely slow. Note also that the fitted 
p, g-distribution is far from perfect which would contribute some amount of error as well. 
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FIG. 11: Left: z = n{q-l) vs logn at K* for the L x L-lattice, L = 6, 8, 10, 12, 14, 16, 32, 64, 
128, 256, 512. The curve is 3 — 1.5 logn — 6 log logn. Right: y = n (p — 1) vs z = n (q — 1) for 
the L X L-lattice, L = 6,8, 10, 12, 14, 16, 32, 64, 128, 256, 512 (512 barely visible near the z-axis). 
The points represent K* . 

The right plot of figure [TT] shows y vs z for a range of temperatures. The points represent- 
ing K* may appear to he on the z-axis but they are are shghtly below it. In the ID-case we 
suspected that there is a limit curve for the high-temperature region, but we suspect that the 
exact data that produced this part of the plot rely on far too small lattices to give any con- 
clusive evidence. Also, the p, g-find algorithm is rather sensitive to noise in this region to be 
useful for sampled data. However, as we said before, this problem goes away once K > K* . 
Figure [12] shows y and z versus K for all the lattices though for the sampled data we only 
show low-temperature data. The red fine is located at K^. = atanh (a/2 — l) ^ 0.44068. 




FIG. 12: y = n{p — 1) vs K (left) and z = n{q — I) vs K (right) for the L x L-lattice, L = 6, 8, 
10, 12, 14, 16, 32, 64, 128, 256, 512. The points represent K* and the red hne is at Kc- The larger 
lattices have their points farther to the right in the plots. 
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C. SD-lattices 



For these lattices we only have exact data for L = 4 and sampled data for L = 
6, 8, 12, 16, 32, 64. The situation is actually somewhat better for 3D-lattices. Figure [T3l 
shows some distributions in the vicinity of K* for L = 32 together with the fitted p, q- 
distributions. For K > K*, just when the distributions become bimodal, the fit is certainly 
less than perfect, but near K* the p, g- approximation is actually rather good. 
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FIG. 13: Magnetisation distributions for the 32 x 32 x 32-lattice (red) and the fitted Pp^g (n, k) 
(blue) Ys k - n/2 for K = 0.2204, K* = 0.22066, K = 0.2210, = 0.2216546 and 0.2220 
(downwards at the y-axis). 

In the left plot of figured!] we show z versus logra at K*. The fitted line through the 
points corresponds to z = —5.3 — logn and is not too bad an approximation. However, 
it was estimated that the growth rate exponent at Kc of the susceptibility is 'y/u = 
1.978 ± 0.009 (assuming 7 = 7' and u = u'). For the magnetisation it was estimated 
fi/v = 0.5147±0.0007. Translated into exponents of n this means 1.657 < 1 +7/3z/ < 1.663 
and 0.8282 < 1 - < 0.8287. If we choose Ai = -5/8 in ffTUSj) and ffTO the first 
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moment exponent would be 53/64 = 0.828125 and 53/32 = 1.65625 for the second moment, 
slightly below the lower bound of the estimate intervals. Choosing Ai = —2/3 would give 
exponents 5/6 = 0.8333 . . . and 5/3 = 1.666 . . . respectively, slightly above the upper bound 
of the intervals. Let us suggest, as an example, that Aq = 6.8, Ai = —2/3, A2 = —6 and 
A3 = in the expression fllOQp . In figure 1X41 the curve use these parameters for z at K* ^ i.e. 
z = 6.8 — (2/3) logn — 6 log logn. Will the points eventually converge to the curve? It would 
take considerably larger lattices to shed any light on this. We also have the problem what 
Aq should be. Using Aq = 6.8 means that the coefficient in f ll09p is about 0.393. Comparing 
the measured a2 with n^/^ gives a factor of roughly 0.16 though the data are certainly far 
from conclusive. Since the distribution fit is not perfect a different constant is perhaps to 
be expected. Also, slow convergence is to be expected here. 



y 




FIG. 14: Left: z = n{q-l) vs logn at K* for the L x L x L-lattice, L = 4,6,8,12,16,32,64. 
The line though the points is —5.3 — logn and the curve is 6.8 — (2/3) logn — 6 log logn. Right: 
y = n (p- 1) vs z for the L X L X L-lattice, L = 4,6, 8, 12, 16, 32, 64 (leftwards) for K > K* . Higher 
values of K when we move downwards left. The red curve is y = 2 u; with w defined by (j59|) . 

The right plot of figure [TH shows y versus z for K > K* . Note the peculiar backwards 
movement of z getting more and more pronounced for larger L. The curves for 16, 32 and 64 
show signs of approaching some limit curve. We don't have data for very low temperatures 
for the smaller lattices though, except for L = 4. The plots in figure [T5l shows y and z versus 
K ioi K > K*. The red lines show location of ~ 0.2216546, found but see also=^ 
for a theoretical estimate of K^. 
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FIG. 15: Left: z = n{q-l)vsK with K > K* for the L x L x L-lattice, L = 4, 6, 8, 12, 16, 32, 64 
(downwards). Right: y = n{p — \) ys K with K > K* for the L x L x L-lattice, L = 
4,6,8,12,16,32,64 (upwards). In both plots the red hne indicates location of Kc and the points 
are the locations of K* . 

D. 4D-Iattices 

In the case of 4- dimensional lattices we have sampled data of magnetisation distributions 
for L = 4, 6, 8, 10, 12, 16. Figure [16] shows some of these magnetisation distributions for 
L = 12 near K* together with fitted p, g-binomial distributions. The fit is quite good, 
considerably better than for 2D and 3D, in the whole range of selected temperatures. Though 
it is hard to distinguish the fitted curves from the magnetisation curves, there is a small 
deviation near the middle. How should z at K* depend on n7 Actually, taking the data at 
face-value they are rather well-fitted to the simple formula z = —6.5 — 0.45 logn. However, 
for the 4D-lattice we have 7 = 7' = l,/3 = l/2 and u = v' = 1/2. This gives that 
l + '-f/du = 3/2 and 1 — P/du = 3/4. Moreover, according tc^^ there should be a correction 
to this. They calculated, using renormalization group techniques, that the susceptibility 
should scale as yTogX near K^. This means that o"2 should scale as n^/^ yTogn. From 
f lll2p we see that we have to choose A2 = —2, with Ai = and A3 = —6, to obtain this. In 
the left plot of figure [T7] we have set z = —1.2 — 2 log logn — 6 log log logn and plotted it 
versus log logn. The curve would then behave as a limit curve rather than as a fitted curve. 
The choice of coefficient Aq = —1.2 is only supported by the human eye as a guide rather 
than any theory and herein lies a problem. With this choice the coefficient of (I112p is about 
0.558. However, dividing the measured a2 at the different K* with n^^"^ A/log n gives values 
close to 0.15. This discrepancy could be due to several sources; e.g. the expression in (I112p 
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FIG. 16: Magnetisation distributions for the 12 x 12 x 12 x 12-lattice (red) and the fitted Pp_g (n, k) 
(blue) vs k- n/2 for K = 0.1490, K* = 0.149255, Kc = 0.149695, K = 0.1500 and K = 0.1505 
(downwards at the y-axis). 

could be incorrect or our data could be suffering from very slow convergence. In the right 
plot of figure [T7| we show y versus z for K > K* together with the curve y = 2w with w 
defined by (1591) . In figure [TS] we show y and z versus K ioi K > K*. The red line is located 
at Kc ^ 0.1496497, estimated in^^. 

E. 5D-lattices 

For the 5-dimensional lattices we have sampled data of magnetisation distributions only 
for L = 4, 6, 8, 10, 12. The distributions in figure[T9]are extremely well fitted by g-binomial 
distributions; it is almost impossible to tell them apart with the naked eye. In five dimensions 
the susceptibility near K^. scales as L^/^, see^. Thus a2 should scale as n^/^ which is exactly 
what we receive when keeping z fixed. So, for z constant we obtain ai oc n^^^ and a2 oc n^/^. 
The left plot of figure [2D] shows z at K* for L = 4, 6, 8, 10, 12. If z approaches a constant 
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FIG. 17: Left: z = n{q-l) vs log log n at K* for the L x L x L x L-lattice, L = 4, 6, 8, 10, 12, 16. 
The curve is —1.2 — 2 log log n — 6 log log log n. Right: y = n(p — 1) vs z = n{q — 1) for the 
L X L X L X L-lattice, L = 4,6, 8, 10, 12, 16 (leftwards). Higher temperatures (low K) begin at the 
upper right part of the plot and with lower temperatures we move down to the left. The red curve 
is y = 2w with w defined by (j59|) . 
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FIG. 18: Left: z = n{q-l)vs K with K > K* ioi the L x L x L x L-lattice, 1 = 4,6, 8, 10, 12, 16 
(downwards). Right: y = n{p — 1) vs K with K > K* for the L x L x L x L-lattice, L = 
4, 6, 8, 10, 12, 16 (upwards). In both plots the red line indicates location of Kc and the points are 
the locations of K* . 

then what is the hmit value? Extracting the hmit z from this plot is futile of course. The 
right plot of figure 120) shows y vs z for the different lattices together with the points K* and 
the curve y = 2w. In figure [21] we show y and z versus K ioi K > K* with an estimated 
marked as a red line. Despite the noise in the plots it seems plausible that z stays essentially 
constant very close to K* (and Kc) and that only y moves. Let us assume this and see where 
this leads us. We employ the moment expressions in section IVIIIl in terms of the parameter a 
to model the behaviour near K*. A normalised first cumulant of the absolute magnetisation 
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FIG. 19: Magnetisation distributions for8x8x8x8x 8-lattice (red) and the fitted Pp^g (n, k) 
(blue) vs A; - n/2 for K = 0.1137, K* = 0.113786, = 0.113914, 0.1143 and 0.1147. 
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FIG. 20: Left: z = n{q-l) vs log n at for the L x L x L x L x L-lattice, L = 4, 6, 8, 10, 12. The 
straight line is constant at z = —9.87. Right: y = n (p—l) vs z = n{q—l) for the LxLxLxLxL- 
lattice, L = 4,6,8,10,12 (leftwards). Higher temperatures (low K) begin at the upper right part 
of the plot and with lower temperatures we move down to the left. The red curve is y = 2w with 
w defined by 

(|M|) I2r?l^ = (Ti/n^/^ should approach gi/ Qq when plotted as a function of a for a fixed z. 
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FIG. 21: Left: z = n{q-l)vsK with K > K* for the L x L x L x L x L-lattice, L = 4, 6, 8, 10, 12 
(downwards). Right: y = n{p — 1) vs K with K > K* for the LxLxLxLx L-lattice, 
L = 4,6, 8, 10, 12 (upwards). In both plots the red Une indicates location of and the points are 
the locations of K* . 

Analogously, the second cumulant (normalised) should behave as 



cr2 - 02 

,3/2 ^ 



(150) 

where the Qm were defined in section IVIIII Note that for a fixed z the Qm now depend only 
on a. The third and fourth cumulants of the absolute magnetisation, divided by respectively 
8n^/^ and 16 n^, quite analogously approach their corresponding limits 
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Through a simple scaling analysis based on our sampled data we have found that the nor- 
malised third cumulant has a limit maximum of about 0.0205 and a minimum of —0.0500. 
The fourth normalised cumulant has a limit maximum of 0.0229 and a minimum of —0.0278, 
based upon our sampled data. Choosing z = —9.87 puts the maximums and minimums of 
the limit curves at appropriate values. Now we identify the coupling K where the minimum 
of the fourth cumulant occurs with the point a where the minimum of the corresponding 
limit curve occurs and likewise for the maximum, thus providing us with a rescaling trans- 
lating K into a. In figure [22] and [23] the first four cumulants are shown together with their 
estimated limit curves for z = —9.87. Indeed the red curve may provide us with a limit. 

Given a lattice size L we denote by K^i^{L) the location of the minimum fourth cumulant 
and by i^Tmaxl-^) the location of the maximum. Analogously for the limit curve, given a 2; we 
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FIG. 22: Normalised first (left) and second (right) cumulants for the LxLxLxLx L-lattice, 
L = 4, 6, 8, 10, 12, 16, 20, 24 (blue) versus a for z = —9.87 together with the limit curve (red). 



FIG. 23: Normalised third (left) and fourth (right) cumulants for the LxLxLxLx L-lattice, 
L = 4,6, 8, 10, 12, 16, 20, 24 (blue) versus a for z = —9.87 together with the limit curve (red). 

denote by a^amiz) and amax(-2) the location of the minimum and maximum fourth cumulant. 
For z = —9.87 we have amin ~ 1.06965 and ctmax ~ 2.51275. A simple scaling projection 
gives that roughly K^^^{L) ^ + and K^^^{L) - K,^UL) ^ OmS/L^/^. Also 

Kc ^ 0.113915, see^^. Thus, in principle at least, the rescaling between a and K is 



However, this kind of expression is somewhat too simplistic to get figure [231 It would 
take higher-order corrections to scaling to produce it but this would probably take a more 
involved numerical study of the 5D-model. Other investigations of the 5D-lattice includes 




K{a) 



(a - amax(2:)) + K^!,x{L) 



(153) 




e.g.ii,^ and2i. 
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XII. CONCLUSIONS 



The magnetisation distribution for the complete graph is exactly described by the p, q- 
binomial distribution, corresponding to the special (or limit) case oi p — q. For balanced 
complete bipartite graphs this is most likely also true in some hmit sense, yet to be made 
precise. Actually, it appears that for most graphs, at least those which are more or less 
regular, the magnetisations are well-fitted by a p, g-binomial distribution for some choice 
of p and q. The exact extent to which the p, g-binomial approximation is good we do not 
yet know (e.g. convergence in moment) nor the exact class of graphs that would satisfy 
this. We have investigated the matter more closely for lattices of dimension one through 
five. In general they are always well-fitted by p, g-binomial distributions for high- and low- 
temperatures but the problems arise near Kc, or rather K* where the distribution changes 
from unimodal to bimodal. 

For the 1-dimensional lattices (having no such bounded K*) the situation is basically al- 
ways that of high temperatures. It seems possible to give expressions for p and q in terms of 
K in this case though we have not done so. For 2-dimensional lattices the distributions near 
K* are least well-fitted by the p, g-binomials but shghtly better fitted in the 3-dimensional 
case. We made theory-based predictions of how z should scale with n near K*. Unfortu- 
nately, scaling is probably very slow, involving logarithms and double logarithms, making it 
near impossible to test the prediction. For 4-dimensional lattices the distributions are clearly 
much better fitted by p, g- binomials, though some discrepancy still remains just above K*. 
For 5-dimensional lattices even this small discrepancy is gone, leaving us perfectly fitted 
(that is, to the human eye) p, g-binomial distributions. In this case the values of z at K* 
should approach a limit value. We estimated this limit and, using this limit value, compared 
the first four normalised cumulants for finite lattices with the (possible) limit curves. 

We described and used a rather simple method to determine p and g given a distribution. 
Possibly this method is not optimal since it simply forces the distribution to be correct at 
a single point rather than providing a good overall- fit. It is also sensitive to noise when 
the distributions are unimodal, thus making it difficult to determine p and g. On the other 
hand it works extremely well for bimodal distributions where the noise sensitivity problem 
vanishes. 

The p, g-binomial coefficients are just a tweaked form of g-binomials, i.e. they are multi- 
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plied by a power of p. It is possible that a different choice of factor would produce better 
results in the case of 2- and 3-dimensional lattices. 

We believe that what is said here for the Ising model also goes for other models, i.e. the 
magnetisation distribution for quantum spin models or for spin-glass models can be modeled 
by p, g-binomial distributions. 
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